A prototype effective-one-body model for non-precessing spinning 
inspiral-merger-ringdown waveforms 



O 

D 

m 



o 
o 

(N 

o 

(N 



X: 



Andrea Taracchini,^ Yi Pan,^ Alessandra Buonanno,^'^ Enrico Barausse,^' ^' * Michael 
Boyle,'' Tony Chu,^ Geoffrey Lovelace,^ Harald P. Pfeiffer,^ and Mark A. Scheel^ 

'Maryland Center for Fundamental Physics & Joint Space-Science Institute, 
Department of Physics, University of Maryland, College Park, MD 20742, USA 
^Radcliffe Institute for Advanced Study, Harvard University, 8 Garden St., Cambridge, MA 02138, USA 
"^Department of Physics, University of Guelph, Cuelph, ON NIG 2W1, Canada 
^Center for Radiophysics and Space Research, Cornell University, Ithaca, NY, 14853, USA 
''Canadian Institute for Theoretical Astrophysics, 60 St. George Street, 
University of Toronto, Toronto, ON M5S 3 118, Canada 
^Theoretical Astrophysics 350-17, California Institute of Technology, Pasadena, CA 91125, USA 

(Dated: February 6, 2012) 

We first use five non-spinning and two mildly spinning (xi — —0.44, +0.44) numerical-relativity 
waveforms of black-hole binaries and calibrate an effective-one-body (EOB) model for non-precessing 
spinning binaries, notably its dynamics and the dominant (2,2) gravitational-wave mode. Then, we 
combine the above results with recent outcomes of small-mass-ratio simulations produced by the 
Teukolsky equation and build a prototype EOB model for detection purposes, which is capable 
of generating inspiral-merger-ringdown waveforms for non-precessing spinning black-hole binaries 
with any mass ratio and individual black-hole spins —1 < X» ^ 0.7. We compare the prototype 
EOB model to two equal-mass highly spinning numerical-relativity waveforms of black holes with 
spins Xi = —0.95, -1-0.97, which were not available at the time the EOB model was calibrated. In 
the case of Advanced LIGO we find that the mismatch between prototype-EOB and numerical- 
relativity waveforms is always smaller than 0.003 for total mass 2O-2OOM0, the mismatch being 
computed by maximizing only over the initial phase and time. To successfully generate merger 
waveforms for individual black-hole spins Xi ^ 0.7, the prototype-EOB model needs to be improved 
by (i) better modeling the plunge dynamics and (ii) including higher-order PN spin terms in the 
gravitational-wave modes and radiation-reaction force. 

PACS numbers: 04.25.D-, 04.25. dg, 04.25.Nx, 04.30.-w 



I. INTRODUCTION 

Coalescing compact-object binary systems (binaries, 
for short) arc among the most promising sources of 
gravitational waves (GWs) for detectors like the U.S. 
Laser Interferometer Gravitational- Wave Observatory 
(LIGO), the British-German GEO, and the French- 
Italian Virgo [1-3]. LIGO and Virgo are undergoing up- 
grades to Advanced configurations [4] , which will improve 
sensitivity by about a factor of 10. A detailed and accu- 
rate understanding of the GWs radiated as the bodies in 
a binary spiral towards each other is crucial not only for 
the initial detection of such sources, but also for max- 
imizing the information that can be obtained from the 
GW signals once they arc observed. 

The matched-filtering technique is the primary data- 
analysis tool used to extract the GW signals from the 
detectors' noise. It requires accurate waveform models 
of the expected GW signals. Analytical templates based 
on the post-Newtonian (PN) approximation [5-8] of the 
Einstein field equations developed over the past thirty 
years accurately describe the inspiraling stage of the bi- 
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nary evolution. In 1999 a new approach to the two-body 
dynamics of compact objects, the so-called effective-one- 
body (EOB) approach, was proposed with the goal of 
extending the analytical templates throughout the last 
stages of inspiral, plunge, merger, and ringdown. The 
EOB approach uses the results of PN theory, black-hole 
perturbation theory, and, more recently, the gravitational 
self-force formalism. It does not, however, use the PN 
results in their original Taylor-expanded form (i.e., as 
polynomials in I'/c), but in a resummed form. 

The EOB formalism was first proposed in Refs. [9, 10] 
and subsequently improved in Refs. [11-13]. Using 
physical intuition and results from black-hole perturba- 
tion theory and the close-limit approximation, Refs. [10, 
13] computed preliminary plunge, merger, and ring- 
down signals of non-spinning and spinning black-hole 
binaries. After breakthroughs in numerical relativity 
(NR) [14-16], the EOB inspiral-merger-ringdown wave- 
forms were improved by calibrating the model to pro- 
gressively more accurate NR simulations, spanning larger 
regions of the parameter space [17-27]. More recently, 
an EOB model for the dominant (2, 2) mode and four 
subdominant modes was built for non-spinning bina- 
ries of comparable-masses [27] and the small mass-ratio 
limit [28]. These results, at the interface between nu- 
merical and analytical relativity, have already had an 
impact in LIGO and Virgo searches. The first searches 



2 



of high-mass and intermediate-mass black-hole binaries 
in LIGO/ Virgo data [29, 30] used the inspiral-merger- 
ringdown templates generated by the EOB model cali- 
brated in Ref. [19], as well as the phenomenological tem- 
plates proposed in Ref. [31]. 

Stellar-mass black holes are expected to carry spins, 
which significantly increases the dimension of the bi- 
nary parameter space. The first EOB Hamiltonian with 
leading-order (1.5PN) spin-orbit and (2PN) spin-spin 
couplings was developed in Ref. [12]. Then, Ref. [13] 
worked out the radiation-reaction force in the EOB equa- 
tions of motion in the presence of spins and computed 
inspiral-merger-ringdown waveforms for generic spinning 
binaries, capturing their main features, including the 
so-called "hang up". Later, Ref. [32] incorporated the 
next-to- leading-order (2.5PN) spin-orbit couplings in the 
EOB Hamiltonian. By construction, in the test-particle 
limit the Hamiltonian of Ref. [32] does not reduce to the 
Hamiltonian of a spinning test particle in the Kerr space- 
time. Moreover, the Hamiltonian of Ref. [32] rewrites 
the EOB radial potential using Fade summation, caus- 
ing spurious poles in some regions of parameter space. 
Nevertheless, the Hamiltonian of Ref. [32] was adopted 
in Ref. [25] to demonstrate the possibility of calibrating 
the EOB model for spinning binaries. 

Since then, substantial progress has been made to- 
wards improving the spin EOB Hamiltonian. Ref. [33] 
worked out the Hamiltonian for a spinning test-particle 
in a generic spacetime, which was used in Ref. [34] to 
derive a spin EOB Hamiltonian having the correct test- 
particle limit. Furthermore, Ref. [34] rewrote the EOB 
radial potential in a way that guarantees the absence of 
poles without employing the Fade summation. As a con- 
sequence, the EOB Hamiltonian of Ref. [34] has desirable 
strong-field circular-orbit features, such as the existence 
of an innermost-stable circular orbit (ISCO), a photon 
circular orbit (or light- ring), and a maximum in the or- 
bital frequency during the plunge. Still preserving these 
properties, the spin EOB Hamiltonian of Ref. [34] was 
recently extended to include the next-to-ncxt-to-leading- 
order (3.5FN) spin-orbit couplings in Ref. [35]. The 
EOB Hamiltonian of Ref. [32] was also recently extended 
through 3.5FN order in the spin-orbit sector in Ref. [36]. 

In the non-conservative sector of the EOB model, the 
radiation-reaction force in the EOB equations of mo- 
tion is built from the GW energy flux, which, in turn, 
is computed from a decomposition of the waveform into 
spherical harmonic {£, m) modes. These modes, instead 
of being used in their Taylor-expanded form, are re- 
summed (or factorized). This factorization was origi- 
nally proposed in Refs. [37, 38] for non-spinning black- 
hole binaries, and was then extended to include spin ef- 
fects in Ref. [39] and higher-order FN spinless terms in 
Refs. [40, 41]. In the test-particle limit, the factorized 
waveforms are known at very high FN order — for ex- 
ample their sum generates the GW energy flux for non- 
spinning binaries through 14FN [41] order and to 4FN 
order in terms involving the black-hole spins. However, 



in the comparable-mass case the GW modes are known 
only at a much lower FN order. Despite the fact that the 
GW energy flux in the comparable-mass case is known 
through 3.5FN [42, 43] and 3FN [44] order in the non- 
spinning and spin-orbit sectors, and 2FN order in the 
spin-spin sector, the GW modes have been computed 
only through 1.5FN order for spin-orbit couplings and 
2FN order for spin-spin couplings [39, 45]. Currently, 
this lack of information in the GW modes is the main 
limitation of our spin EOB model, and, as we will see, it 
affects the performance of the model for prograde orbits 
and large spin values. 

In this paper, we build upon the past success in an- 
alytically modeling inspiral-merger-ringdown waveforms 
through the EOB formalism, and develop a prototype 
EOB model for non-precessing spinning black-hole bina- 
ries that covers a large region of the parameter space and 
can be used for detection purposes and future calibra- 
tions. More specifically, we adopt the EOB Hamiltonian 
derived in Refs. [34, 35], the GW energy flux and factor- 
ized waveforms derived in Refs. [38, 39], and calibrate the 
EOB (2.2) dominant mode to seven NR waveforms: five 
non-spinning waveforms with mass ratios 1, 1/2, 1/3, 1/4 
and 1/6 [27] and two equal- mass non-precessing spinning 
waveforms of spin magnitudes 0.44 [46]. We combine the 
above results with recent small-mass-ratio results pro- 
duced by the Teukolsky equation [28] and build a proto- 
type EOB model for inspiral-merger-ringdown waveforms 
for non-prcccssing spinning black-hole binaries with any 
mass ratio and individual black- hole spins — 1 < Xi ^ O-'''- 
For Xi ^ 0.7, although the EOB dynamics can be evolved 
until the end of the plunge, the EOB (2,2) mode peaks 
too early in the evolution, where the motion is still qua- 
sicircular. As a consequence, we cannot correct the EOB 
(2,2) mode to agree with the NR (2,2) mode peak using 
non-quasicircular amplitude coefficients. This limitation, 
which also affects the small- mass-ratio limit results [28], 
is caused by the poor knowledge of FN spin effects in the 
GW modes and makes the prototype EOB waveforms 
unreliable for Xi ^ 0.7. Two NR waveforms with nearly 
extremal spin magnitudes [47, 48] became available to 
us when we were finishing calibration of the spin EOB 
model. We use them to examine the limitations of the 
spin prototype EOB model, and extract from them useful 
information for future work. 

The paper is organized as follows. In Sec. II, we de- 
scribe the spin EOB model used in this work, its dy- 
namics, waveforms, and adjustable parameters. Sec- 
tion III A discusses the alignment procedure used to com- 
pare EOB and NR waveforms at low frequency, and the 
statistics used to quantify the differences between the 
waveforms. We then calibrate the EOB model to the 
NR waveforms in Sec. IIIB. In Sec. IV, we combine the 
results of Sec. Ill A with those of Ref. [28] to build a 
prototype EOB model that interpolates between the cal- 
ibrated EOB waveforms and extends them to a larger 
region of the parameter space. We also investigate how 
this prototype EOB model performs with respect to two 
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NR waveforms with nearly extremal spin, which were not 
used in the calibration. Finally, Sec. V summarizes our 
main conclusions. In Appendix A we explicitly write the 
factorized waveforms used in this work, including spin 
effects. 



II. EFFECTIVE-ONE-BODY DYNAMICS AND 
WAVEFORMS IN THE PRESENCE OF SPIN 
EFFECTS 

In this section, we define the spin BOB model that we 
will later calibrate using NR waveforms. Henceforth, we 
use geometric units G ~ c = 1. 

In the spin BOB model [12, 32, 34-36] the dynamics of 
two black holes of masses mi and m2 and spins S'l and 
S^2 is mapped into the dynamics of an effective particle 
of mass fj. = mi 7712/ {mi + m2) and spin moving in 
a deformed Kerr metric with mass M = mi + 7712 and 
spin S'Korr- The position and momentum vectors of the 
effective particle are described by R and P, respectively. 
Here, for convenience, we use the reduced variables 



r = 



R 



P 

p= — 



(1) 



Since we will restrict the discussion to spins aligned 
or anti-aligned with the orbital angular momentum, we 
define the (dimcnsionlcss) spin variables Xi = 
Xi rnf L, where L is the unit vector along the direc- 
tion of the orbital angular momentum. We also write 

S'Kerr = XKerr-M^L. 



A. The efFective-one-body dynamics 

In this paper we adopt the spin BOB Hamiltonian pro- 
posed in Refs. [33-35]. The real (or BOB) Hamiltonian 
is related to the effective Hamiltonian through the 
relation 



i/rcal = M^froal = MWl + 2l/ 



off 



- 1 



M, 



(2) 



where Hoff describes the conservative dynamics of an ef- 
fective spinning particle of mass ^ and spin S* moving 
in a deformed Kerr spacetime of mass M and spin S'Korr- 
The symmetric mass-ratio v = fi/M acts as the deforma- 
tion parameter. Through 3.5PN order in the spin-orbit 
coupling, the mapping between the effective and real spin 
variables reads [34, 35] 



^Kcrr — Si - 

S* = ^ 
mi 



Si 



mi 
m2 



.(1) 



(3a) 
(3b) 



where A^V and a'-^} are the 2.5PN and 3.5PN spin-orbit 
terms given explicitly in Bqs. (51) and (52) of Ref. [35]. 
They depend on the dynamical variables r and p, the spin 



variables Si, and on several gauge parameters. These pa- 
rameters are present because of the large class of canon- 
ical transformations that can map between the real and 
effective descriptions. Their physical effects would cancel 
out if the PN dynamics were known at arbitrarily high 
orders; since this is clearly not the case, the gauge param- 
eters can have a noticeable effect [35] and may in princi- 
ple be used as spin BOB adjustable parameters. In this 
paper however, we set all gauge parameters to zero and 
introduce a spin BOB adjustable parameter at 4.5PN or- 
der in the spin-orbit sector by adding the following term 
to Bq. (3b) 



.(3) 



dso V fm2 



mi 



Si 



mi 
m2 



(4) 



Here d^o is the spin-orbit BOB adjustable parameter. 
The effective Hamiltonian reads [34] 



13' p, + aJl + Y^p^PJ + Qiip) + + 

//V MM 

(5) 

where the first two terms are the Hamiltonian of a non- 
spinning test particle in the deformed Kerr spacetime, 
a, and 7*^ arc the lapse, shift and 3-dimcnsional 
metric of the effective geometry and Qi{p) is a non- 
geodesic term quartic in the linear momentum intro- 
duced in Ref. [49]. The quantities Hso and i?ss in 
Bq. (5) contain respectively spin-orbit and spin-spin cou- 
plings that are linear in the effective particle's spin S* , 
while the term -~l/{2Mr^){r'^5^^ - 3r'r^)S*S* is the 
leading-order coupling of the particle's spin to itself, 
with S''^ being the Kronecker delta. More explicitly, us- 
ing Ref. [34] we can obtain Hso and H^s by inserting 
Bqs. (5.31), (5.32), Bqs. (5.47a)-(5.47h), and Bqs. (5.48)- 
(5.52) into Bqs. (4.18) and (4.19); a, /3* and 7*^ are given 
by inserting Bqs. (5.36a)-(5.36c), Bqs. (5.38)-(5.40) and 
Bqs. (5.71)-(5.76) into Bqs. (5.44)-(5.46). We will eluci- 
date our choice of the quartic term Q4{p) at the end of 
this section, when introducing the tortoise variables. 

Following Ref. [25], we introduce another spin BOB 
adjustable parameter in the spin-spin sector. Thus, we 
add to Bq. (5) the following 3PN term 



mi 



Si 



mi 
m2 



(Si + S2 



(6) 



with dgg the spin-spin BOB adjustable parameter. For 
what concerns the non-spinning BOB sector, we adopt 
the following choice for the BOB potentials At and 
entering a, (3i and 7^ (sec Bq. (5.36) in Ref. [34]). The 
potential At is given through 3PN order by 

At(u) = ^A„(it), (7a) 

A„(w) = A(M)+xLrrW', (7b) 

A{u) =^ l-2u + 2uu'^ + v (—-—-kA ^7^(7c) 

\ 3 32 / 
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where u = 1/r. Reference [34] suggested rewriting the 
quantity Au{u) as 

Au{u) = A„(ii) [l + Ao + log (l + Ai M + A2U^ 
+A3u^ + Aiu^)] , (8) 

where A^ with i = 1,2,3,4 are exphcitly given in 
Eqs. (5.77)-(5.81) in Ref. [34], and 



sionless form by [25] 




1 



„EOB 
'H,- 



(9a) 



(9b) 



Here, are radii reducing to those of the Kerr event 

and Cauchy horizons when the EOB adjustable param- 
eter K goes to zero. The logarithm in Eq. (8) was in- 
troduced in Ref. [34] to quench the divergence of the 
powers of u at small radii. Its presence also allows the 
existence of an ISCO, a photon circular orbit (or light- 
ring), and a maximum in the orbital frequency during 
the plunge. The reason for modeling A„(u) with Eq. (8) 
instead of using the Pade summation of A„(u), as pro- 
posed in Ref. [32], is threefold. First, we did not want 
to use the Fade summation of A„(m) because Ref. [25] 
found that for certain regions of the parameter space 
spurious poles can appear. Second, although wc could 
have applied the Fade summation only to A{u) and used 
the Fade potential A{u) calibrated to non-spinning wave- 
forms in Ref. [27], we want to take advantage of the good 
properties of the potential (8) during the late inspiral, as 
found in Ref. [34]. Third, we find it useful to develop a 
variant of the EOB potential so that in the future we can 
test how two different EOB potentials (both calibrated 
to NR waveforms at high frequency) compare at low fre- 
quency. 

Furthermore, for the potential A^ at 3FN order enter- 
ing the EOB metric components (5.36) in Ref. [34], we 
choose 

Ar{u) = At{u)D-'^{u) , (10a) 
D^\u) = l + log[l + 6i/u^ + 2(26-3i/)i^u^]. 

(10b) 

Once expanded in FN orders, the EOB Hamiltonian 
(2) with the effective Hamiltonian defined in Eq. (5) 
and the spin mapping defined in Eqs. (3a) and (3b), 
reproduces all known FN orders — at 3FN, 3.5FN and 
2FN order in the non-spinning, spin-orbit and spin-spin 
sectors, respectively — except for the spin-spin terms at 
3FN and 4FN order, which have been recently com- 
puted in Refs. [50-57]. Furthermore, in the test-particle 
limit the real Hamiltonian contains the correct spin-orbit 
couplings linear in the test-particle spin, at all FN or- 
ders [33, 34]. 

Let t = t/M. In terms of the reduced Hamiltonian 
i/ioai, the EOB Hamilton equations are given in dimen- 



dr - 9i?rcal 
-J = {r, iJroalj = , 

dt op 



(11a) 



5 = {p,^real}+^=-^+i^, (lib) 

dt or 

where ^ denotes the non-conservative force that ac- 
counts for radiation- reaction effects. Following Ref. [13], 
we use ^ 



dE 



i^fl\r X p\ dt 



(12) 



where ft = M\r x f|/r^ is the dimensionless orbital fre- 
quency and dE /dt is the GW energy flux for quasicircular 
orbits obtained by summing over the modes (^, m) as 



dE Cl^ \ ^ X ^ 

dt Stt ^ ^ 

1=2 m=0 



(13) 



Here TZ is the distance to the source, and simply elimi- 
nates the dominant behavior of /i£m- We sum over pos- 
itive m modes only since |/i£,m| = |/i£,-m|- Expressions 
for the modes are given in the next section. In this 
paper, we restrict the calibration to non-precessing bi- 
naries, and thus we omit the Hamilton equations of the 
spin variables. 

It was demonstrated in previous work [37, 58] that by 
replacing the radial component of the linear momentum 
Pr = (jJ ■ r)/r with p^*, which is the conjugate momen- 
tum of the EOB tortoise radial coordinate r* , one can 
improve the numerical stability of the EOB equations 
of motion. This happens because Pr diverges when ap- 
proaching r^'^^ while Pr* does not. In this paper we fol- 
low the definition of the EOB tortoise radial coordinate 
in Appendix A of Ref. [25].^ However, when applying the 
tortoise coordinate transformation to the quartic term in 
Eq. (5), we get [25] 



Q4(p)(x4?^(^'+XkcJ', 



(14) 



which clearly diverges &i r — r^^^ . As in the non- 
spinning case [27, 37, 58], we neglect contributions higher 
than 3FN order and rewrite Eq. (14) as 



Q4(p)(x^(r2+xLj\ 



(15) 



which is well behaved throughout the EOB orbital evo- 
lution. 



^ The over-dot stands for d/dt. 
Note that all the formulas in Appendix A of Ref. [25] are written 
in physical dynamical variables, namely R and P, while here we 
use reduced variables r and p. 
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Lastly, we integrate the EOB Hamilton equations. In 
order to get rid of any residual eccentricity when the EOB 
orbital frequency is close to the initial frequency of the 
NR run, we start the EOB evolution at large separation, 
say 50M, and use the quasispherical initial conditions 
developed in Ref. [13]. We stop the integration when the 
orbital frequency 51 reaches a maximum. 



B. The efFective-one-body waveforms 

Following Refs. [24-27, 37] we write the inspiral-plunge 
modes as 



7 insp-plungc 



im 1 



(16) 



where the hj^ are the factorized modes developed in 
Refs. [37-39], and the Nim are non-quasicircular (NQC) 
corrections that model deviations from quasicircular mo- 
tion, which is assumed when deriving the /if,^ . The fac- 
torized modes read 



{N,, 



I A(e) rp 



(17) 



where e is the parity of the waveform. All the factors 
entering the /if^ can be explicitly found in Appendix A. 
We emphasize here again that despite the fact that the 
GW energy flux in the comparable-mass case is known 
through 3PN order in the spin-orbit sector [44] , the spin- 
orbit couplings in the factorized (or PN-expanded) modes 
have been computed only through 1.5PN order [39, 45]. 
This limitation will degrade the performances of our spin 
EOB model for prograde orbits and large spin values, as 
already observed in the test-particle limit in Refs. [28, 39]. 
To improve the knowledge of spin effects in the GW 
modes, Refs. [25, 39] added spin couplings in the test- 
particle limit through 4PN order in the factorized wave- 
forms. However, since the mapping between the Kerr 
spin parameter in the test-particle limit and the black- 
hole spins in the comparable-mass case is not yet unam- 
biguously determined [34, 35], and since we do not have 
many NR spinning waveforms at our disposal to test the 
mapping, wc decide not to include here the spinning tcst- 
particle-limit couplings in the factorized waveforms com- 
puted in Ref. [39]. Wc have checked before performing 
any calibration that EOB models with or without test- 
particle spin effects (with Kerr spin parameter XKorr) give 
similar performances. 

In all the calibrations of the non-spinning EOB model, 
two EOB adjustable parameters were needed to calibrate 
the EOB Hamilton equations — for example Refs. [26, 27] 
used the 4PN and 5PN order coefficients in the EOB po- 
tential A{r). As discussed in the previous section, for 
the EOB model adopted in this paper, the EOB non- 
spinning conservative dynamics depend so far only on the 
adjustable parameter K. We introduce a second EOB ad- 
justable parameter in the non-conservative non-spinning 
EOB sector by adding a 4PN order non-spinning term 
in P22 and denote the coefficient of this unknown PN 



term by pgt'' [^cc Eq. (A8a)]. This adjustable parameter 
enters the EOB Hamilton equations through the energy 
flux defined in Eq. (13). 

As shown in Ref. [27], the NQC corrections of modes 
with {£,771) (2,2) have marginal effects on the dy- 
namics. Also, our goal in this work is to calibrate only 
the (2, 2) mode, so in the following we set Ni„i = 1 for 
(£,m) ^ (2,2). We have^ 




X exp 




(18) 



where af^^ (with i = 1...5) are the (real) NQC amplitude 
coefficients and (with i — 1...4) are the (real) NQC 
phase coefficients. We will explain in detail how these 
coefficients are determined at the end of this section. 

The EOB merger-ringdown waveform is built as a lin- 
ear superposition of the quasinormal modes (QNMs) of 
the final Kerr black hole [10, 17, 19, 21, 23, 24, 59], as 



N-l 



I mcrgcr-RD 



(t) = A22n e 



-i(722n [i — in 



(19) 



ri=0 



where N is the number of overtones, A22n is the complex 
amplitude of the n-th overtone, and a22n = ^22n — i/T22n 
is the complex frequency of this overtone with positive 
(real) frequency uj22n and decay time T22n- The com- 
plex QNM frequencies are known functions of the mass 
and spin of the final Kerr black hole. Their numerical 
values can be found in Ref. [60]. The mass and spin 
of the final black hole, Mf and a/, can be computed 
through analytical phenomenological formulas reproduc- 
ing the NR predictions. Here, we adopt the formulas 
given in Eq. (8) of Ref. [61] and in Eqs. (1) and (3) 
of Ref. [62]. We notice that the formula for the final 
mass in Ref. [61] was obtained using numerical simula- 
tions of small-spin black-hole binaries with mildly un- 
equal masses. As a consequence, the formula is not very 
accurate for the large-spin, unequal-mass binaries con- 
sidered in this paper. However, other formulas available 
in the literature are either very accurate but only valid 
for equal-mass binaries [63], or have not been yet exten- 
sively tested against NR simulations [64, 65]. Thus, for 
the time being we stick with Eq. (8) of Ref. [61], but we 
plan to construct a better formula in the future using all 
recent data in the literature. 

Furthermore, we follow the hybrid matching procedure 
of Ref. [27] to fix the N complex amplitude coefficients 



^ Note that in Ref. [28] the Ni^^ were written in terms of physical 
dynamical variables, rather than the reduced variables used here. 
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A22n in Eq. (19). We set up N complex linear equa- 
tions by imposing that the inspiral-plunge and merger- 

1 r 7 inspiral-plunffc j 7 merger- RD 

rmgdown waveiorms n.22 ^ ^ " and /i22 co- 
incide on iV — 2 points (evenly sampled over a range 
[^match - ^^match. ^matchl) ^nd that thdr time derivatives 

^^inspiral-plunge /^mcrge.-RD ..i^.i^C at tH,^,,^ At^,^^ 

and i^atcii- previous works, we introduce the EOB 

adjustable parameter At'^^^^^ which describes the size of 
the comb over which we impose continuous and smooth 
matching in order to determine the ringdown waveform. 

In Refs. [24, 27, 28], pseudo QNMs (pQNMs) were 
proposed and applied to moderate the rise of the EOB 
GW frequency during the merger-ringdown transition — 
for example Sec. IIC of Ref. [27] discussed in some de- 
tail the advantage of using pQNMs for higher-order GW 
modes. In this paper, we find it useful to introduce a 
pQNM for the (2,2) mode. Therefore, we choose N = 8 
in Eq. (19) and replace the highest overtone in the sum- 
mation with this pQNM. 

Finally, we build the full inspiral-plunge-merger- 
ringdown EOB waveform by joining the inspiral-plunge 
waveform ft,2T''"'''"'''""^''(0 ^^d the merger-ringdown 
waveform /i22°'^''''^^(i) at the matching time i^atcii 

+ /i™-^°(t)(?(t-iL\tch)- 

In Fig. 1, we summarize how the inspiral-plunge- 
merger-ringdown EOB waveform is constructed. Beyond 
the ISCO, the quasi-circular inspiral waveform is followed 
by a short plunge waveform ^ where the radial motion is 
no longer negligible and NQC corrections quickly become 
important. The plunge ends roughly when the effective 
particle in the EOB description crosses the light-ring, 
which, in the non-spinning case, coincides approximately 
with the peak of EOB orbital frequency fi and waveform 
amplitude \h22\- Until this moment, the GW radiation in 
the EOB description is obtained directly from the motion 
of the effective particle. After this moment that we iden- 
tify as the merger, the direct emission of GWs from the 
effective particle is strongly attenuated and filtered by 
the potential barrier formed around the newborn black 
hole. Thus, in the EOB description the merger-ringdown 
waveform is no longer obtained from the motion of the 
effective particle, but it is built through a superposi- 
tion of QNMs. This procedure of constructing the full 
EOB waveform, in particular replacing the direct emis- 
sion with a superposition of QNMs beyond the light ring, 
was first proposed in Refs. [10, 13] for non-spinning and 
spinning comparable-mass black-hole binaries. It was in- 
spired by the close limit approximation [66] and results 
in Refs. [67, 68] where it was observed that once the 



* The number of gravitational-wave cycles during the plunge scales 
roughly as [10]. 




50 100 150 200 



FIG. 1: We show in the spacetime diagram {t,r*) the trajec- 
tory of the effective particle in the EOB description (black 
solid line in the left part of the diagram) and the EOB 
(2,2) gravitational mode (red solid oscillating line) for an 
equal-mass non-spinning black-hole binary. Although we only 
need to evolve the EOB trajectory until the orbital frequency 
reaches its maximum ("light ring"), the model's dynamics al- 
lows the trajectory to continue to negative r* (short-dashed 
black line in the left part of the diagram). The blue dashed 
lines represent t ± r* = const, surfaces and ingoing/outgoing 
null rays. The EOB (2,2) mode is a function of the retarded 
time t — r* , plotted here orthogonal to t — r* = const, sur- 
faces, at a finite t + r* distance. The two outgoing null rays 
are drawn at the t — r* retarded times when the EOB particle 
crosses the EOB ISCO and light-ring radii, respectively. The 
shaded green area is a rough sketch of the potential barrier 
around the newborn black hole. 



radially infalling particle is inside the potential barrier 
which peaks around the light ring, the direct gravita- 
tional radiation from the particle is strongly filtered by 
the potential barrier. Part of the energy produced in the 
strong merger-burst remains stored in the resonant cav- 
ity of the geometry, i.e., inside the potential barrier, and 
what is released outside is just the ringdown signal. The 
non-linear scattering of GW radiation (tails) against the 
curvature potential of the newborn black hole also con- 
tributes to the merger-ringdown waveform. Currently, 
in the EOB description the merger-ringdown waveform 
is effectively the tail of a ^-function impulse at merger. 
When spin effects are present, the overall picture de- 
picted in Fig. 1 survives, but with some differences due 
to the fact that the EOB light-ring position, peak of the 
orbital frequency ft and waveform amplitude I/122I can 
be displaced in time [28]. We notice that the physical 
picture of the merger-ringdown that emerged from the 
studies in Refs. [66-68] and was incorporated in the EOB 
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description in Rcfs. [10, 13], has also recently motivated 
the hybrid approach of Refs. [69, 70]. 

We now continue our detailed review of how the EOB 
waveform is built and discuss how we fix the NQC coeffi- 
cients in Eq. (18). Since we do not expect spin effects in 
the NQC correction until 1.5PN order in either amplitude 
or phase, the coefficients a^*^^ with i = 1,2 and b'l'' with 
i = 1,2 only depend on v, while af^'^ with i = 4, 5 and 
b^^^ with z = 3, 4 are functions of v linearly proportional 
to the spins Xi,2- The coefficient Og^^ is given by the sum 
of a non-spinning term (dependent only on v) and a spin- 
ning term (proportional to the spins Xi,2)- In Sec. Ill we 
first calibrate the non-spinning waveforms, and then the 
spinning ones. Thus, we determine the ten coefficients in 
Eq. (18) in two steps. First, we set Xi — X2 ^ 0, thus 
= (with i = 4, 5) and b'l^^ = (with i = 3, 4) 
and calculate the values of the five NQC coefficients a^^^ 
(with i = 1,2, 3) and b'l^^ (with i = 1, 2) by imposing the 
following five conditions [27, 28]: 

1 . Let ipcak be the time at which the EOB orbital fre- 
quency reaches its peak. Then, the peak of the 
EOB (2, 2) mode must happen at the matching 



time imatch 



peak 



Ay-22 



dt 



that is 



0, 



(21) 



where Atp^^j. is an EOB adjustable parameter, 
which will be specified in Sec. III. We note that 
in Ref. [28] the quantity At^^^^y. was computed by 
comparing the times at which the Teukolsky (2,2) 
mode and the EOB orbital frequency reach their 
peaks. This was possible because the EOB trajec- 
tory was used in the Teukolsky equation to evolve 
the dynamics. However, in the NR simulation, we 
do not know what Atp^^^j, is, because the EOB dy- 
namics does not determine the NR dynamics. 

2. The amplitudes of the NR and EOB (2, 2) modes 
are the same. 



|,EOB^.n , A. 



22 N 
peak J 



~ l''22 v'-peak 



)l 



(22) 



3. 



The curvatures of the amplitudes of the NR and 
EOB (2, 2) modes are the same. 



dt'^ 



(23) 



4. The GW frequencies of the NR and EOB (2,2) 
modes are the same. 



"^22 l^^pcak "T ^^-pcak) — ^22 V-pcak) 



(24) 



The time derivatives of the GW frequency of the 
NR and EOB (2, 2) modes are the same. 



dt 



peak po 



dt 



(25) 



We summarize in Table I all the NR-input values that we 
use in the right-hand side of Eqs. (22)-(25). After the 
five non-spinning NQC coefficients have been computed, 
we plug them back into the EOB dynamics through the 
energy flux, start a new EOB evolution, generate a new 
EOB (2, 2) mode, and calculate new NQC coefficients. 
We repeat this procedure until the values of the NQC 
coefficients converge. Then, when calibrating spinning 
waveforms, we set a^*^^ and fe^*^^ (with i = 1,2), as well 
as the non-spinning part of flg^^, to the values just cal- 
culated for Xi = X2 = 0, and apply the five conditions 
above in an iterative way, obtaining the final coefficients 
a'l^^ (with « = 3, 4, 5) and (with i = 3, 4). Note that 
in order to generate GW templates, this procedure can 
be computationally expensive, since to generate one EOB 
(2, 2) mode one has to evolve the dynamics a few times. 
The current computational cost of generating an EOB 
waveform long enough for the LIGO bandwidth varies 
between a fraction of a second to a few seconds,^ depend- 
ing on the masses. The iterative procedure can increase 
this cost by a factor of a few. 

In order for the NQC coefficients to be effective in cor- 
recting the EOB mode peak, the latter has to occur in 
a region where the radial motion is comparable to or at 
least ^30% of the tangential motion. Such a condition 
is in principle not a necessary requirement for the EOB 
model to work. In fact, the radial motion is expected to 
be strongly suppressed for almost extremal black holes, 
at least in the test-particle limit, since the ISCO coincides 
with the horizon for X ~ 1 [71]- However, if the factor- 
ized (2,2) mode, given by Eq. (17), differs substantially 
from the NR (2,2) mode because of the lack of high-order 
spin-orbit terms, the inability of the NQC coefficient to 
change the waveform during the plunge at high spins may 
prevent the EOB model to work properly. This is because 
the NQC coefficients cannot artificially compensate the 
missing higher-order spin orbit terms in the waveforms, 
as they partially do at low spins. In fact, we will see 
that this problem arises for Xi }^ 0.7, making the EOB 
prototype waveforms unreliable for large positive spins. 

We list in Table II all the EOB adjustable parameters 
that we exploit in this work to calibrate the EOB model 
to NR simulations. 



III. EFFECTIVE-ONE-BODY CALIBRATION 

In this section, we calibrate the EOB model using 
seven NR waveforms, namely five non-spinning wave- 
forms of mass ratios q = TO2/mi = 1,1/2,1/3,1/4 
and 1/6 and two equal- mass spinning waveforms with 
Xi = X2 = -H0.43655 and xi = X2 = -0.43757. The 
calibration is achieved by minimizing the amplitude and 



The time is measured by running a code that is not optimized in 
speed on a single CPU. 
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TABLE I: Exact NR-input values used in the right-hand side of Eqs. (21)-(25) to cahbrate the EOB inspiral-plunge waveforms. 



Xl = X2 


1 




1/2 



1/3 



1/4 



1/6 



1 

+0.43655 


1 

-0.43757 


l"-22,peakl 


0.3940 


0.3446 


0.2855 


0.2403 


0.1810 


0.3942 


0.3935 




-10.3 


-8.8 


-6.9 


-5.5 


-3.9 


-7.7 


-12.4 




0.3593 


0.3467 


0.3324 


0.3218 


0.3084 


0.3989 


0.3342 




11.3 


10.5 


9.6 


8.9 


8.1 


11.2 


10.7 



phase differences between the NR and EOB (2, 2) modes 
over the six EOB adjustable parameters: K, dgo and 
dss in the EOB conservative dynamics, and pjt^ ^^peak' 

in the EOB waveforms (see 



PQNM , 

Table II). 



pQNM 
^"22 



A. Alignment of EOB and NR waveform 

When calibrating NR and EOB waveforms, we first 
align the waveforms at low frequency following the pro- 
cedure of Refs. [24, 25, 27]. This procedure consists of 
minimizing the square of the difference between the NR 
and EOB (2,2)-mode phases (j)22^ and (/)f°^, integrated 
over the time window (^1,^2)1 



if2°^i + M 



/>?2^(0] rft. 



(26) 



with respect to the time shift to a-nd phase shift , where 
it is understood that ^fl^^ is computed for a chosen set of 
adjustable parameters. The time window (^1,^2) should: 
(i) begin as early as possible, where the NR and EOB 
GW-phase evolutions agree best, (ii) begin late enough 
to avoid the junk radiation present in the numerical sim- 
ulation, (iii) be long enough to average over numerical 
noise, and (iv) extend from peak to peak (or trough to 
trough) over an integer number of oscillations in the GW 
frequency, which are caused by the residual eccentricity 
in the numerical initial conditions. In Table III, we list 
our choices of (^1,^2) f^or the seven numerical waveforms 
at our disposal. Each time window extends through 10 
eccentricity oscillation cycles in the numerical frequency 
evolution. 

Let (/)o and to be the alignment parameters. Then, 
we define the phase and relative amplitude differences 
between the EOB and NR (2,2) modes as follows: 



/EOB 
"22 



and 



\ A J 



- 1 



(27) 



(28) 



We then define the global phase and relative amplitude 
differences over a time window (tijts) with 



A^giobai = max \A(j}{t)\, 

t6(tl,t3) 



(29) 



TABLE IL Summary of adjustable parameters of the spin 
EOB model considered in this paper. The values of the 
EOB adjustable parameters used in this paper are given in 
Eqs. (33), (34), (35), (36), (37), and (39). In addition, the 
NQC parameters a^^^ and 6^^^ are fixed from NR-input val- 
ues through Eqs. (21)~(25). 



EOB dynamics 
adjustable parameters 



EOB waveform 
adjustable parameters 



K 

dso, rfss 



A .22 A ,22 

"''match 1 "'■peak 
pQNM pQNM 



and 



AA 



= max 

global *e(ti,i3) 



/AA 



(30) 



In the following, when measuring the difference between 
NR and EOB inspiral-plunge waveforms we set = 
^match ' while whcn we measure the difference between full 
inspiral-merger-ringdown waveforms we use ^3 = tend, 
where tend is chosen as late as possible into the ringdown 
stage, but before numerical errors due to gauge effects 
become noticeable [24]. We list the values of tfj^atch ^^^^ 
tend for the seven NR waveforms in Table III. 



B. Procedure to calibrate the EOB adjustable 
parameters 

Recently, Ref. [28] computed the waveforms in the 
small mass-ratio limit by evolving a time-domain Teukol- 
sky equation in which the source term is evaluated using 
an EOB trajectory. It was found that there exists a time 
difference between the Teukolsky (2, 2)-mode amplitude 
peak and the EOB orbital-frequency peak. This differ- 
ence is parametrized by the quantity At^^^^^ introduced 
in Eq. (21). Table III in Ref. [28] lists this difference as a 
function of the Kerr spin parameter: for non-spinning 
and retrograde cases — 3M < Atp^j^j^. < 1.6Af, while for 
prograde cases At^^^^ decreases quickly as function of 
the spin. Let us consider XKcrr, which explicitly reads 



XKo 



(l-2i/) XS + V1-41/XA, 



(31) 
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and also define 



X = Xs + XA ■ 



\-2v 



(32) 



where xs,A = (xi ±X2)/2. For an equal-mass, equal-spin 
binary {v = 1/4, xi = X, X2 = x) we have XKcn- = x/2, 
while in the test-particle limit we have XKgit = X (that is 
the spin parameter of the background spacetime) . There- 
fore, inspired by the results in the test-particle limit, we 
assume here that for an equal-mass, equal-spin binary 
^^pcak depends on the black-hole spins through x- Ex- 
plicitly we choose 



A ,22 
^'peak 



-2.5M- 1.77A/I 



if X < , 
if X > , 



(33) 



which models qualitatively Table III in Ref . [28] . Follow- 
ing Refs. [24, 25, 27], we calibrate the EOB adjustable pa- 
rameters in two steps. These steps are performed for each 
of our seven calibration NR waveforms separately, result- 
ing in seven sets of calibration parameters. First, for each 
of the NR waveform at our disposal, we use At'^^^^ in 
Eq. (33), insert the NR- input values from Table I into 
Eqs. (21)-(25), solve them iteratively for the NQC coef- 
ficients, and calibrate K, pi^-i (or c'so and dss if spins are 
present) by minimizing Eq. (29) with = t 



match' 



This 

process provides us with the EOB inspiral-plunge wave- 
form. Second, to obtain the EOB merger-ringdown wave- 
form, we calibrate the size of the comb At^J^^^j^ and the 
pQNM (complex) frequency by applying Eq. (29) with 
^3 = tend- As in Ref. [27], we find that a constant value 
for the comb size, notably 

At^ltch = 7.5Af, (34) 

gives a very good performance for all the different mass 
ratios and spins. A detailed study of the pQNM (com- 
plex) frequency has revealed that the best result is ob- 



tained when w??^^ 



■^22 

uj^^^M/Mf at 



match 



lies between the GW frequency 
and the frequency of the least- 



damped QNM L1J220, and when r. 



pQNM 
22 



is (not much) 



shorter than T22o- Specifically, we use the simple choice 



UJ. 



pQNM 
22 



,EOB 



{t 



22 \ 
match/ 



M 



- ^220 



dQNM 
'22 



Y^r220: 



(35a) 
(35b) 



for all different mass ratios and spins. Before ending 
this section, we discuss in more detail how we carry out 
the calibration of the parameters K, ^22^, for the non- 
spinning sector, and the parameters dso, (isS; for the 
spinning sector. 



TABLE III: We list the parameters t\, t2 entering the align- 
ment procedure defined in Eq. (26), and the parameter 
(both thatch s-iid tend) entering the computation of waveforms' 
differences in Eqs. (29) and (30). 



Xl = X2 


1 




1/2 



1/3 



1/4 



1/6 



1 

+0.43655 


1 

-0.43757 


ti/M 

^match /^^ 
tcnd/M 


820 
2250 
3943 
3990 


770 
2255 
3729 
3770 


570 
1985 
3515 
3560 


670 
1985 
3326 
3370 


870 
2310 
4892 
4940 


800 
2150 
3367 
3410 


610 
1850 
2402 
2430 



K{iy) and i'^) by a global minimization of A(/)giobai 
and (AA/A)giobai [as defined in Eqs. (29) and (30) us- 
ing = i^atch] with respect to the unknown coefficients 
entering the Kliy) and P22 {'^) polynomials. However, as 
in previous studies [26, 27], we find a strong degeneracy 
among the EOB adjustable parameters, when calibrating 
each mass ratio separately. The degeneracy is partially 
broken when we combine all the available mass ratios 
together, but it is not completely lifted. In particular, 
different choices of K{h') and p^J (v) lead to EOB mod- 
els that can match equally well with NR waveforms. We 
are thus relieved from a rigorous yet expensive global 
search and follow a simplified procedure to find satis- 
factory K{v) and p^22{^)- First, we locate two points 
(0.8154, -35) and (1.188, -20) in the K-p'§ plane where 
A^giobai < 0.1 rad and (AA/^)giobai < 0.1 for g = 1 and 
q = l/Q {v ^ 0.25 and v = 0.1224), respectively. We 

then determine a linear function P22' (^) by imposing that 
p^2^(0.25) = -35 and p^2^ (0.1224) = -20, leading to 



-5.6 - 117.61/. 



(36) 



At (7 = 1/2,1/3 and 1/4, we choose p'^22 according to 
Eq. (36) and determine the value of K that minimizes 
A^giobai and a range of K values that satisfy A(/)giobai < 
0.1 rad. 

We now have a complete set of calibration parameters 
for each of our non-spinning NR waveforms. In order to 
obtain calibration parameters that interpolate between 
the NR waveforms, we build a least-squares fit quadratic 
in v against these K values. By construction, we fix 
two of the three free parameters in the fit by requiring 
that in the test-particle limit K{v) reproduces the ISCO 
shift of Refs. [34, 72, 73] and that the optimal equal- 
mass value i^(0.25) is recovered exactly. Even with these 
two constraints and just one free parameter to fit, the 
residuals are within 1% (see Fig. 2). We find 



Calibrating non-spmmng waveforms 



K{v) = 1.447 - 1.715 V - 3.246 



(37) 



In general, the adjustable parameters K and ^22"* de- 
pend on the mass ratio and we assume that they are poly- 
nomial functions of v. In principle, we should determine 



Finally, since the iterative procedure to compute the 
NQC coefficients through Eqs. (21)-(25) can be expen- 
sive, we have parametrized them through quadratic fits. 
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0.05 0.1 0.15 0.2 0.25 

V 

FIG. 2: We show the quadratic fit in v for the adjustable 
parameter K. This parameter is cahbrated using the five non- 
spinning NR waveforms, assuming fi~^^{y) in Eq. (36). The 
error bars are determined by the intersection of the contours 

of A^giobai = 0.1 rads with p^^-iiy) for each mass ratio con- 
sidered. 



finding rather small residuals. Explicitly, we obtain 

a^^^ = -12.68 + 75.421^- 106.6 (38a) 

02"" = 101.5 - 757.3 J/ -f 1473 (38b) 

03"" = -107.7 + 857.6 1/- 1776 (38c) 

h\^^ = -1.464+ 12.82 1/- 60.10 1/^ (38d) 

62'' = 7.477 - 85.26;/ + 353.3 (38e) 

2. Calibrating spinning waveforms 

When calibrating the EOB inspiral-plunge waveforms 
to the two NR equal-mass, equal-spin waveforms at our 
disposal (xi = X2 = +0.43655 and Xi = X2 = —0.43757), 
we use the non-spinning EOB adjustable parameters K 
and Eqs. (37)-(36), and calibrate the spinning 

EOB adjustable parameters dso and dss- We reach this 
goal by building contour plots in the plane dso-dss for 
AtZ-giobai in Eq. (29) with = t^^tch- We find that the 
contours of A</)giobai = 0.2 rads associated with the two 
NR spinning waveforms intersect each other for the fol- 
lowing choice of the adjustable parameters 

dso = -69.5, dss = 2.75. (39) 

Note that when computing the spinning NQC coeffi- 
cients, we use the NQC coefficients parametrized in 
Eq. (38), and solve iteratively the five conditions (21)- 



(25) for a'l^^ {i = 3, 4, 5) and b'l^^ {i = 3, 4).** 

IV. A PROTOTYPE EFFECTIVE-ONE-BODY 
MODEL FOR NON-PRECESSING SPINNING 
WAVEFORMS 

We now build on the results of Sec. Ill, and also 
on recent outcomes of small mass-ratio simulations pro- 
duced by the Teukolsky equation [28] , to construct a self- 
contained set of prescriptions to generate EOB inspiral- 
merger-ringdown waveforms in a larger region of the pa- 
rameter space (^, X1JX2) of the binary. 

A. Interpolating the EOB model outside the 
domain of calibration 

Since we only have seven NR waveforms at our disposal 
(and just two of them with spins), when extending the 
EOB model to regions of the parameter space without 
NR waveforms, we are forced to make assumptions on 
the behavior of the adjustable parameter Atp^^^j^ and the 
NR-input values in Table I. In this work we assume that 
the 3 dimensional space (i^, xi, X2) can be treated as the 2 
dimensional space {1/, x)- [Note that v & [0, 1/4] and x € 
[—1,1].] More specifically, given a binary described by 
the parameters (^, Xi!X2) having in general xi 7^ X2j we 
consider an auxiliary equal-spin binary with parameters 
(j/, X,x)i where x is defined as in Eq. (32). With this 
choice, the auxiliary binary has the same value of XKcn- as 
the original binary. We stress that the auxiliary binary is 
used only to extend the EOB adjustable parameters and 
the NR-input values to regions of the parameter space in 
which we do not have NR results. Of course the EOB 
dynamics and waveforms are computed for the original 
binary, not the auxiliary one. 

Thus, in the prototype EOB model, the EOB ad- 
justable parameter Atp^g^j^ in Eq. (33) is evaluated us- 
ing for X the value from Eq. (32). To compute the 
spinning NQC coefficients in the prototype model, we 
need to prescribe the input values in the right-hand side 
of Eqs. (22)-(25) using the parameters of the auxiliary 
binary. We proceed as follows. We only have knowl- 
edge of the NR-input values at merger for a few re- 
gions of the {ly, x) parameter space. We can obtain the 
NR-input values along the curve {i^ = 0, x) from the 
Teukolsky waveforms of Ref. [28]. In particular, both 
l^22!poakl and 9t2|/i^2^peakl are set to (since they are 
proportional to ^), while for w^j^pg^^j, and i^2?peak we 
use the data in Table V of Ref. [28]. We can extract 
the peak information along the curve {v = 1/4, x) from 
the three equal-mass waveforms used in the calibration 



Note that the NQC coefRcient a^^^ is solved for twice, first in 
the non-spinning calibration and then in the spinning one. 
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TABLE IV: Fits of the NR-input values /^^ that are used to build the global fits in 
Eq. (42) for the test-particle and equal-mass limits. 

f^^' Curve Fit 



iu = 0,x) 
22.p«akl (j,^x/4,x) 0.3961 



r2a2ii.NR 



IVl Of |'l22,pcak 



(i^ = l/4,x) lO-" X (-1.007 + 0.5415x) 

^NR {i^ = 0, x) 0.2758 - 0.08898 log(l - x) 

'^'^'^'"'^ {u = 1 /4, x) 0.3604 + 0.08242X + 0.02794x^ 



■ NK 

J« '^22, peak 



{u = 0, x) 10"^ X [5.953 + (0.7199 + 1.210x) log(l - x)] 

(i^ = l/4,x) 0.01113 



of this paper, together with the two nearly extremal spin 
cases xi = X2 = -0.94905 and xi = X2 = +0.9695 
(not used for the calibration of the adjustable parame- 
ters dso c?ss)j which we will discuss in Sec. IV D. 
Along the curve {ly, x = 0) we can use the NR-input val- 
ues of the non-spinning waveforms from Refs. [27, 28]. 
In Table IV we list the fits for each NR-input value 

G mi.J^d?\hl%.J^^^t^^^^ in the 

test-particle and equal mass limits. Along the non- 
spinning profile, fits quadratic in i/ give a good descrip- 
tion of the exact NR-input values, hence we assume that 
the dependence of on i' is quadratic as well and has 
the simple form 

/N«(i., x) = C2(x) + ci (x) >^ + co(x) . (40) 

We can fix two of the coefficients Ci by imposing that 
the test-particle limit and equal-mass cases are exactly 
recovered when v ^ and v = 1/4, respectively. We 
can fit the third coefficient to the exact NR-input values 
along the non-spinning direction. This means that the 
fits along the non-spinning profile are not exactly recov- 
ered by the global fits /^^(i', x), but we find that the 
residuals are negligible. Explicitly, we fit ci in the fol- 
lowing expression 

/"«(;., 0;ci) = {16[,ri^(l/4,0)-/Ni^(0,0)]-4ci}z.2 
+ ciz. + ,ri^(0,0), (41) 

and denote the fitted value with ci. Finally, we extend 
the result outside the non-spinning profile assuming that 
the global fit reads 

/^« (i/, x) = {16[/N«(l/4, x) - /^« (0, x)] - 4ci} 

+ cii^ + /NR(0,x). (42) 

In Table V we list the values of ci for the four NR-input 
values that are needed to compute the right-hand sides 
in Eqs. (22)-(25). 

Having in hand ^t^^^ak the NR-input values, we 
complete the construction of the prototype EOB model 
by fixing the EOB adjustable parameters K, P22^ ^-^id 



TABLE V: Fitted values of ci for the four NR-input values 
as defined in Eq. (42). 





I^NR 

l"-22,peak 


M^^ll/l^^peakl 


A T, ,NR 

iV-J 0^22, peak 


jt2,- NR 

"^22, peak 


Cl 


1.355 


-2.5 X 10"-' 


0.1935 


0.01204 



rfso, ds,s to the values in Eqs. (37)- (36) and (39), re- 
spectively, employing the pQNM (complex) frequency in 
Eq. (35), the comb size in Eq. (34), and the NQC coeffi- 
cients in Eqs. (38). 

To test the robustness of the construction of the quan- 
tity /^'^(j^, x): we study how the spinning NQC coeffi- 
cients change across the plane (;/, x)- We focus on bina- 
ries with xi = X2 = X- We compute iteratively the NQC 
amplitude coefficients a^*^^ (with i = 3,4,5) for different 
mass-ratios in the range 1/100 < 9 < 1 and for different 
spins in the range — 1 < Xi ^ 0.7 {i = 1,2). Typically, 
we get convergence of the NQC coefficients within five it- 
erations. Unfortunately, we cannot span larger, positive 
values of Xi since the NQC corrections tend to diverge 
as the spin magnitude grows in the prograde case. The 
reason is that they become less effective in reshaping the 
EOB (2,2) peak as prescribed by the fits f^^{v, x)- This 
happens because the peak of the EOB (2,2) mode oc- 
curs too early in the evolution when the orbital motion 
is still quasicircular. Hence the NQC coefficients must be 
very large to compensate for the small values oi p^^ / {r^) 
and be able to reshape the EOB (2,2) amplitude around 
the peak in a satisfactory way. As discussed earlier, this 
would not be a problem in principle if higher-order spin- 
orbit terms were known in the factorized waveforms, but, 
as a result of the lack of knowledge of those, our EOB 
prototype waveforms are reliable only up to Xi ^ 0.7. 



B. Performance for non-spinning waveforms 

In Figs. 3 and 4 we show how the inspiral-merger- 
ringdown EOB waveforms computed according to the 
prescriptions of Sec. IV A compare with the NR wave- 
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FIG. 3: Comparison of the NR and BOB (2, 2) mode for q=l, 
Xi = X2 = 0. In the upper panels we show the comparison 
between the real part of the two waveforms, zooming into the 
merger region in the upper right plot. In the lower panels we 
show the dephasing and relative amplitude difference over the 
same time ranges as the upper panels. A vertical dashed line 
marks the position of the NR amplitude peak. The dotted 
curves are the NR errors. 




1000 2000 3000 4000 4800 4850 4900 495' 




0.4 
0.3 
0.2 
0.1 
0.0 
-0.1 
-0.2 



A(|) ; 
AA/A 

















1000 2000 3000 4000 4800 4850 4900 495 
(t - r J/M (t - r J/M 



0.4 
0.3 
0.2 
0.1 
0.0 
-0.1 

r 



FIG. 4: Same as in Fig. 3 but for g = 1/6, xi = X2 = 0. 



forms for two representative mass ratios q = 1,1/6. In 
general, for all the non-spinning vifaveforms we find that 
the dephasing is typieally within 0.1 rads up until if,Jatch 
(merger time) and always within 0.2 rads when ineluding 
the ringdown stage. The figures also show in dotted lines 
the NR phase and amplitude errors obtained by combin- 
ing the extrapolation and resolution errors in quadrature. 
We notiee that the EOB and NR amplitudes' agreement 
is remarkably good up to the merger time, while during 
the ringdown the relative amplitude difference may grow 
up to about 15%, approaching the estimated NR error. 



In Ref. [27] the authors calibrated a different version 
of the non-spinning EOB model to the same set of non- 
spinning NR waveforms used in this paper, the main dif- 
ference between the two EOB models being the choice of 
the EOB potential A{r), as we discussed in Sec. II A. We 
find that the difference between the EOB inspiral-merger- 
ringdown waveforms and the NR waveforms in Ref. [27] 
is comparable to and for some mass ratios marginally 
worse than what we have achieved in this work using 
the prototype EOB model. The only noticeable qualita- 
tive difference is that the phase error of the prototype 
EOB model accumulates more slowly during the merger- 
ringdown transition because of the introduction of the 
pQNM in the (2, 2) mode. We point out that the in- 
clusion of the pQNM (complex) frequency in the EOB 
merger-ringdown waveform is not strictly needed for the 
non-spinning case, but we use it even in this case for 
uniformity with the spinning sector, where the pQNM 
frequency is instead crucial. 

We can quantify the differences between NR and EOB 
waveforms by computing the mismatch (MM), as defined 
in Eq. (43) of Ref. [27], which is one minus the overlap be- 
tween two waveforms, weighted by the noise spectral den- 
sity of the detector and maximized over the initial time, 
phase and binary parameters. If we use an Advanced 
LIGO noise curve, named ZERO_DET_HIGH_P in Ref. [4], 
wc obtain that the MM, maximizing only over the initial 
phase and time, is always smaller than 0.001 when the 
binary total mass varies between 20Mq and 2OOM0. For 
these total masses, the NR waveforms start in band. We 
taper them using the Planck-taper window function [74] 
to reduce numerical artifacts. The width of the window 
function is set to the length of NR waveforms, ranging 
from O.35(M/2OM0) to 0.65{M/20Mq) seconds. The 
window function smoothly rises from to 1 in the first 
0.0625 seconds and falls from 1 to in the last 0.0125 sec- 
onds. We restrict the MM integration to the frequency 
band for which NR waveform is available. 



C. Performance for spinning waveforms 

In Figs. 5 and 6 we present the results of the pro- 
totype EOB model for the two spinning waveforms at 
our disposal. We observe that the choice (39) gives a 
larger dephasing for Xi = X2 = +0.43655 than for xi = 
X2 = —0.43757 or the non-spinning runs. In fact at the 
merger time the dephasing for the Xi = X2 = +0.43655 
waveform grows beyond the NR error. For the ampli- 
tude, we instead get a similar performance, on the same 
level as the other runs. The worse performance of the 
Xi = X2 = +0.43655 waveform can be explained in the 
more relativistic nature of this run. In fact, in this case 
the EOB ISCO moves to smaller radial separations as 
the spin parameter x increases towards positive values 
(aligned runs) . On the other hand, for negative values of 
X (anti-aligned runs) the EOB ISCO moves outwards to a 
less relativistic regime and one expects a better behavior 
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of the EOB model. This expectation is confirmed by the 
cahbration of the Xi = X2 = —0.43757 run, for which we 
find that very good performances can be achieved in large 
regions of the EOB adjustable parameter space. Fig. 5 
shows that in this case the dephasing is well within the 
NR error at the merger time. For these spinning wave- 
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FIG. 5: Same as in Fig. 3 but for g = 1, Xi = X2 = -0.43655. 
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FIG. 6: Same as in Fig. 3 but for g = 1, Xi = X2 = -fO.43756. 



forms, we obtain that the MM, maximizing only over the 
initial phase and time, is always smaller than 0.003 when 
the binary total mass varies between 20 Mq and 200Mq. 



D. Performance for nearly extremal spin 
waveforms 



Here we compare the EOB waveforms of the prototype 
model developed in Sec. IV A, against two equal-mass 



NR waveforms with nearly extremal spins: Xi = X2 = 
-0.94905 and xi = X2 = +0.9695 [47, 48]. We stress that 
these NR waveforms were not used when calibrating the 
spin EOB adjustable parameters o?so f^nd dgg in Eq. (39). 
The only information that we used from these two nearly 
extremal spin waveforms was their NR-input values when 
building the fits /^^(j^, x)- 

As already discussed, when the spins are anti-aligned, 
the EOB ISCO moves towards larger radial separations, 
so that the binary is less relativistic throughout its or- 
bital evolution as compared to the aligned configura- 
tions. Therefore, we expect that in this case the EOB 
model is more effective. The results in Fig. 7 for the case 
Xi = X2 = —0.94905 confirm this expectation. The de- 
phasing grows up to about 2 rads during the ringdown, 
while the relative amplitude difference grows up to about 
40%. Despite the large phase difference at merger, we 
find that, even without maximizing over the binary pa- 
rameters but only the initial phase and time, the MM 
is always smaller than 0.005 for systems with total mass 
between 2OM0 and 2OOM0. 

For the case Xi = X2 = +0.9695, which is outside the 
domain of validity of our prototype EOB model, we can- 
not successfully run the NQC iterations, since the NQC 
corrections are so large that they cause a divergent se- 
quence of NQC coefficients. Nonetheless, we deem it 
interesting to generate the EOB inspiral-plunge wave- 
form where only the non-spinning NQC coefficients a^*^^ 
(i = 1, 2, 3) and b'l'^'^ [i = 1, 2) are used and compare it to 
the NR waveform. In Fig. 8 we show how our waveform 
performs. We notice that the NR waveform is very long, 
almost 50 GW cycles. The phase difference between the 
EOB and NR waveforms is smaller than 0.04 rads over 
the first 20 GW cycles, and then grows up to 0.18 rads 
during the subsequent 10 GW cycles and it becomes 0.9 
rads when 10 GW cycles are left before merger. The 
fractional amplitude difference is only 3% when 10 GW 
cycles are left before merger. 

It is worth emphasizing that although our prototype 
model is not yet able to generate merger-ringdown wave- 
forms for spins larger than +0.7, nevertheless, as the com- 
parison with the nearly extremal case Xi = X2 = +0.9695 
has proven, the Hamiltonian of Refs. [34, 35] and the 
resummed flux of Refs. [38, 39] can evolve the EOB dy- 
namics in this highly relativistic case beyond the orbital- 
frequency's peak, until r « 1.9M, without encountering 
unphysical features. This suggests that relevant strong- 
field effects are well grasped by the EOB dynamics and 
waveforms [34, 35, 38, 39], at least as far as the NR runs 
used in this paper are concerned. Moreover, the large 
amplitude difference causing the NQC iteration to break 
down for large, positive spins was already observed in 
Refs. [28, 39] where it was pointed out that it is im- 
portant to improve the modeling of spin effects in the 
EOB waveform amplitude. Finally, as observed above, 
the breaking down of the NQC procedure in this highly 
relativistic case, although not a problem in principle if 
higher-order spin-orbit terms were known in the factor- 
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izcd waveforms, is due to the fact that the peak of the 
EOB (2,2) mode occurs too early in the orbital evolution 
where non-quasicircular orbit effects are still negligible. 
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FIG. 7: Same as in Fig. 3 but for g = 1, Xi = X2 = -0.94905. 
This NR waveform was not used to calibrate the adjustable 
parameters dso and dss- Alignment between the NR and 
EOB waveforms was performed using Eq. (26), with ti — 
860 M and t2 = 2470 M. 
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FIG. 8: Same as in Fig. 3 but for g = 1, Xi = Xa = H-0.9695 
and only the inspiral portion. This NR waveform was not used 
to calibrate the adjustable parameters dso and dss- Also, 
in the aligned case our prototype EOB model only covers 
Xi,2 J; 0.7. Note that in this plot we do not include spinning 
NQC corrections in our EOB waveform. Alignment between 
the NR and EOB waveforms was performed using Eq. (26), 
with ti = 1170 M and ti = 2790 M. 



V. CONCLUSIONS 

Using the EOB spin Hamiltonian in Refs. [34, 35], the 
factorized waveforms in Refs. [38, 39], and the adjustable 
parameters in Table II, we have developed a prototype 
EOB model for non-precessing spinning black-hole bina- 
ries that can be used for detection purposes in LIGO and 
Virgo searches and employed for future calibrations [75]. 
The prototype model is built by first calibrating the 
EOB adjustable parameters against five non-spinning 
waveforms with mass ratios 9=1, 1/2, 1/3, 1/4, 1/6 and 
two equal-mass, equal-spin NR waveforms with moderate 
spins. Then, those results, at the interface with NR, are 
combined with recent results at the interface with black- 
hole perturbation theory [28]. The resulting prototype 
EOB model interpolates between calibrated points in the 
binary parameter space, and generates inspiral-merger- 
ringdown waveforms with any mass ratio and individual 
spin magnitudes — 1 < Xi ^ 0.7. This EOB model has 
been implemented in the freely available LIGO Algorithm 
Library (LAL) [76] with the model name "SEOBNRvl".^ 

Wc found that the EOB waveforms generated with the 
prototype model agree with the NR waveforms used to 
calibrate them within '-^O.l rads at merger for the non- 
spinning sector, and within ~0.15 rads at merger for 
the spinning sector. In terms of amplitude differences 
at merger, both non-spinning and spinning runs agree 
to within 5%. The MMs for Advanced LIGO computed 
by maximizing only with respect to the initial phase and 
time are always smaller than 0.003 for binaries with total 
masses between 20 Af© and 2OOA/0. 

Wc also compared the prototype EOB model to two 
equal-mass, equal-spin NR waveforms of black holes with 
nearly extremal spins, notably Xi — —0.94905, -1-0.9695. 
Those NR waveforms were not part of the original set of 
waveforms used to calibrate the EOB model. Wc found 
that for the anti-aligned case the prototype EOB model 
performs quite well for detection purposes, with MMs 
smaller than 0.003 without maximizing over the binary 
parameters, but only on initial phase and time. In the 
aligned case, which is highly relativistic due to a spin 
as large as -1-0.9695 (outside the range of vahdity of our 
prototype model) , we compared the inspiral-plunge wave- 
form for 40 GW cycles and found a dephasing of ^0.8 rad. 
During the last 10 GW cycles before merger the dephas- 
ing grows up to several radians. This non-satisfactory 
performance during plunge and merger for large, posi- 
tive spins is not surprising. In our prototype spin EOB 
model the factorized modes [39] used in the radiation- 
reaction force generate spin couplings in the GW energy 
flux at a PN order much lower than what is known to- 
day. In fact, the GW energy flux is currently known 



Two non-spinning EOB models arc also available in LAL, "EOB- 
NRvl" and "EOBNRv2" , which were calibrated to NR wave- 
forms in Refs. [19, 27]. 
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through 3PN order in the spin-orbit sector^ [44] and 2PN 
order in the spin-spin sector. However, the —2 spin- 
weighted spherical harmonics that are used to build the 
factorized waveforms employed in this paper are known 
only through 1.5PN order in the spin-orbit sector [45]. 
Moreover, the performance wc found for large spin values 
and prograde orbits confirms what was already found in 
Ref . [28] , where EOB waveforms in the test-particle limit 
could be calibrated to Teukolsky-type waveforms only up 
to a Kerr spin value of ~ -I- 0.7. For larger spin values, 
the factorized waveforms start deviating from the exact 
ones even before reaching the IS CO [28, 39]. 

The prototype spin EOB model can be improved in the 
future in different directions. First, the choice of the spin 
EOB adjustable parameters done in Sec. II was rather ar- 
bitrary and assumed that all gauge parameters that enter 
the spin EOB conservative dynamics are zero. Of course, 
it would have been difficult to carry out a more sophisti- 
cated study in this work considering that we had at our 
disposal only two equal-mass, equal-spin NR waveforms. 
When several more spin NR waveforms will be avail- 
able, the spin EOB parameters (together with the non- 
spinning ones) should be explored and calibrated simul- 
taneously against all the available NR waveforms. Sec- 
ond, it is urgent to compute higher-order PN spin-orbit 
terms in the -2 spin weighted spherical harmonics and in 
the factorized modes, thus making the EOB spin model 
reliable also for large, positive spins, i.e., for Xi > 0.7. 
Third, the spin EOB Hamiltonian at 3.5PN order used 
in this paper predicts for large, positive spins that the 
position of the peak of the EOB orbital-frequency varies 
non-monotonically as function of the spin and lies in a re- 
gion which is not very relativistic. It would be important 
to correct this behavior calibrating the gauge parameters 
present in the spin EOB Hamiltonian. Fourth, recent re- 
sults in Refs. [73, 77-79] at the interface between PN the- 
ory and the self-force formalism, have allowed Ref. [80] to 
compute the non-spinning EOB potentials at all orders 
in PN theory, but linear in the symmetric mass ratio v. 
These new results will be incorporated in the future to 
improve the non-spinning conservative dynamics of the 
prototype EOB model, and will be extended to include 
spin effects. 
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Appendix A: Explicit expressions of the factorized 
modes 



Using results from Refs. [24-27], we write here the ex- 
plicit expressions of the factorized modes employed in 
Sec. II B. Even though we calibrated only the (2,2) mode, 
we will provide expressions for all the modes up to £ = 8, 
because they enter the computation of the energy flux in 
Eq. (13). 

The terms h^,^''^ in Eq. (17) are the Newtonian modes. 



They read 



e'^-f^"tlc...(.)v^ir— (|,$), (Ai) 

where TZ is the distance from the source; the F^'" (0, $) 
are the scalar spherical harmonics; we use = v^'^ 
with 
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where p$ = ]r x p). The functions n^^^^^ and C(+t{v) in 
Eq. (Al) read 
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and 



2 2 



+ (-i)^+M'i + ix/T^ 



The function S^\}i^ in Eq. (17) is an effective source term 
that in the circular-motion limit contains a pole at the 
EOB light ring. It is given in terms of the EOB dynamics 
as 



(A4) 



Reference [48] found that the tail spin-orbit terms in the energy 
flux at 3PN order dominate all the other spin-orbit contributions 
and improve the agreement with NR waveforms. 



^'=V,Pr.,P$), e = 0, 
^ off = I'o , e = 1 ; 



(A5) 
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where vq = Cl^/^. The factor Tf„ in Eq. (17) resums the 
leading order logarithms of tail effects, it reads 
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where ro = 2A//^ [39]. 
In what follows we define 



5m \ 



mi — 7712 



Also we use eulerlog„(7;Q) = 'Je + log 2 + logm 
(A7a) l/21og7;Q, with 'Je being the Euler constant. 
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where is a non-spinning EOB adjustable parameter, which is determined through the calibration of the non- 
spinning NR waveforms, 
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where the 1.5PN spin terms in Eqs. (A9a) and (A9c) correct typos in Eqs. (42e) and (42e) of Ref. 
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176 1/2 - 850 Z/ + 375 3 410 833 4 
390 (2 J/ - 1) ~ 709 800 ' 

21 980 z/'"'- 104 930^2^84 679^- 15 828 2 7187914 4 

13 650 (51/2 - 5^+1) ~ 15 526 875 ' 

8 !/2 - 626 ;/ + 319 2 31877 4 
390 (2 1/ - 1) ~ 304 200^^^ ' 



15 526 875^ 



(Alia) 
(AUb) 
(Allc) 
(AUd) 
(Alle) 



P66 



P65 



P6i 



P63 



P62 



Pel 



273 z/'^ - 861 ^2 + 602 ^ - 106 2 

84(5^2 _ 5 J, ^ ;L) " 

220i/3-910i/2 + 838t/-185 2 
144 (3z/2-4i/+l) ' 

133j/3-581t/2 + 462^-86 
84(5^2 _ 5 1) ^ 

156 z/3- 750 ^2 ^ 742 J,. 



1025 435 
659 736 



476887 
659 736 



169 



144 (3 1/2 -4 J/ +1) 
49 1/3 - 413 :/2 + 378 1/ - 74 2 

84(5j/2- 5i/+ 1) 
124 ^3 _ 670 ^2 + 694 ^ - 161 

144 (3 ^2 _ 4 4^1) 



817991 
3 298 680 



/Q ' 



(A12a) 
(A12b) 
(A12c) 
(A12d) 
(A12c) 
(A12f) 



13801/3 - 4963t/2 + 4246i/ - 906 , /» o n 

P77 = 1 H : — ^1 ^ Va , (A13a 

^' 714(3i/2 -4i/+ 1) ^ ^ 

. 6104J/4 - 29351^3.^ 37828^2 „;Lg^85j^_^ 2144 2 

= ^ + 1666(7.3-14.2 + 7.-1) ' 

804.3 - 3523.2 + 3382.- 762 , , 
- 1 + 714 (3.2- 4. + 1) ' 



Pf4 = l 



P73 = 1 



P72 = l 



P71 



- 217959J/3 + 298872i/2 _ i31805j/ + 17756 2 
14994 (7i/3-14t/2 + 7i^-l) ^ 
420J/3 - 25631/2 + 2806j/ - 666 , 

714 (3i/2 _4j, + i) 
32760:/4 - 190239^/3 _^ 273924z/2 _ 123489^/ + 16832 

14994 (7j/3-14z/2_^7i/-1) 
228J/3 - 20831^2 _^ 25181/ - 618 2 



714 (3i/2-4i/ + l) 
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(A13d) 
(Al3e) 
(A13f) 
(A13g) 



P^7 



P86 = 



PS5 



Pal 



P83 



PS2 



Psi 



3482 - 26778i/ + 646591/^ - 534451/^ + 12243;/^ 2 

2736(-l + 7t/- 14z/2 + 7i/3) ' 
23478 - 1540991/ + 309498i/2 - 2075501/^ + 38920:/^ 

18240(-1 + 61/ - lOz/2 + 4i/3) 
1002 - 7498t/ + 17269t/2 - 13055j/3 + 2653;/^ 2 
912(-l + 7i/- 14:/2 + 7z/3) 

4350 - 280551/ + 54642^/2 - 345981/^ + 6056;/^ 

3648(-l + 6i/-10i/2 + 4t/3) 
2666 - 194341/ + 42627i/2 - 289651/^ + 4899:/^ 

2736(-l + 7z/ - 14z/2 + 7z/3) 
20598 - 131059i/ + 249018i/2 - 149950j/3 + 24520:/"^ 

18240(-1 + 6i/ - 10j/2 + 4i/3) 
2462 - 175981/ + 37119i/2 - 228451/^ + 3063i/^ 

2736(-1 + 7j/-14i/2 + 7i/3) 
20022 - 1264511/ + 236922t/2 - 138430^/^ + 21640:/^ 



18240(-1 + 61/ - 10j/2 + 4i/3) 
Finally, we give the explicit expression of the phase term 



4 



7 

S22 = -^(^ Hicul 



V 315 
I 



-24^/^. 



(Al4a) 
(A14b) 
(Al4c) 
(A14d) 
(Al4e) 
(Al4f) 
(A14g) 
(A14h) 

(A15) 
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